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Abstract. We present universal continuous variable quantum compu- 
tation (CVQC) in the micromaser. With a brief history as motivation 
we present the background theory and define universal CVQC. We then 
show how to generate a set of operations in the micromaser which can 
be used to achieve universal CVQC. It then follows that the micromaser 
is a potential architecture for CVQC but our proof is easily adaptable to 
other potential physical systems. 



1 Introduction 

Analogue computation has a long-running history, from the invention of the 
Astrolabe [1] for plotting the heavens in around 200 BC, through the slide rule 
and mechanical differential analyser to more modern electronic devices. Analogue 
computation is less well developed than its digital counterpart, but offers many 
opportunities both in the theoretical advancement and physical realisation of 
computers [2-4]. 

Quantum mechanics is a more recent invention, being conceived and devel- 
oped since the early 20th century. Much of the original research in quantum 
mechanics used continuous variable systems, such as operations on the positions 
and momenta of particles. It would thus appear that quantum mechanics offers 
a breeding ground for new theories of continuous variable (CV) computation. 
However, perhaps inspired by the prevalence of classical digital computation, 
most of the research into quantum computation is aimed at discrete variables, 
in the form of qubits [5-7] . 

Lloyd and Braunstcin [8] laid the groundwork for continuous variable quan- 
tum computation (CVQC) in 1999. Almost immediately research went into im- 
plementing CV algorithms, such as analogues of Grover's, Deutsch & Jozsa's 
and Shor's algorithms [9-11], and investigating the general structure of compu- 
tation and simulation with CVs. Some of these looked at implementing discrete 
computation embedded in CV systems [12, 13] and others looked at using the 
physics of CV systems directly to implement CV computing [14, 15]. 

There are two major schemes for CVQC: encoding the information in infinitely- 
squeezed states, such as the position eigenstates; or encoding the information in 
Gaussian states, such as the quantum coherent states of light. Previous work 



favours the first approach but we argue that the second is more appropriate for 
a realistic view and implementation of CVQC. 

In this paper we will show how to achieve universal CVQC in a particular 
experiment, the micromascr. Computations and universality are described in 
Section 2 with recipes for universality in two different encodings of variables. The 
micromaser itself is described in Section 3 along with the background physics. 
Our results are laid out in Section 4, which is how to achieve universal CVQC 
in the micromaser and in Section 5 we give our plans for further work. 



2 Universal CVQC 

The notion of universality is important in all branches of computation theory. We 
call our system a universal computer for our purposes if it can perform any com- 
putation for which wc wish to use it. In discrete variable quantum computation 
(DVQC), universal computation is being able to achieve any unitary operation 
on the states encoding the variables. Since realising any unitary operation on a 
continuous variable would require an infinite number of parameters, for continu- 
ous variable quantum computation (CVQC) we restrict ourselves to exponentials 
of Hermitian polynomials on the space of continuous variables. This is sensible as 
unitary operations are usually considered to be the result of applying a Hamil- 
tonian for a period of time and Hamiltonians are Hermitian polynomials ^ . 
We encode the CV information in the eigenstates of some continuous-spectrum 
operator and computations are embodied as physical manipulations which cor- 
respond to operations on the eigenstates. This is the definition of continuous 
variable quantum computation that we employ. To encode our variables and to 
describe the physical modes involved we make use of the quadrature operators 
X and p, which are orthogonal in the sense that [x,p] = i, up to a real normali- 
sation constant. Any polynomial in x and p can be generated given a certain set 
of available operators as stated by Lloyd & Braunstein [8]: 

"Simple linear operations on continuous variables, together with a non- 
linear operation and any interaction suffices to enact to an arbitrary 
degree of accuracy arbitrary Hermitian polynomials of the set of contin- 
uous variables." 

In terms of operators, the simple linear operations are {±^, ±p}, a non-linear 
operation might be the Kerr Hamiltonian Hkbtt = {x^ + P^) and an interaction 
is to couple modes together, e.g. two mode squeezing or the sum gate described 
later in this section. The way in which polynomials arc generated is as follows: 
given Hamiltonians A and B, for some small time 5t 
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^ A polynomial / in the position and momentum operators is said to be Hermitian 
iff it is its own adjoint: P = f 



and 
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So we can generate the Hamiltonians ±i A, ±B and ±A ± B from and B 
very easily to arbitrary fidelity. Note that we are converting between the space of 
unitary operators and the space of Hamiltonians to generate our Hamiltonians 
since unitary evolution is how nature evolves. A non-linear operation C is one of 
order cnibc'd or higher in x and p, meaning that when it is commuted with another 
operator H, an operator of higher order in x and p than H is obtained. Non- 
linear operations can't be efficiently simulated on a classical discrete variable 
computer [14]. Recursively then, polynomials of any order can be generated. 

While this gives us universality by our definition of obtaining any Hermitian 
polynomial in x and p, there is some choice left in how to encode the CV quanti- 
ties into the states. The set of continuous variable states is generally considered 
to belong to one of two distinct classes. One consists of variables encoded into 
a set of infinitely squeezed states - eigenstates of a quadrature operator - and 
the other is the set of Gaussian states. The infinitely squeezed states are a limit 
case of the Gaussian states but the two have slightly differing sets of universal 
operations. The two are not equivalent but lead to the same set of computable 
functions since we are looking for any Hermitian polynomial in the operators. 

We choose to encode our information in finitely squeezed Gaussian states 
because infinitely squeezed states are somewhat unphysical. To achieve universal 
computation we need a non-linear operation (order cubed or higher) and the 
following list of linear operations (in order: displacement; fourier transform; 1- 
mode squeezing; 2-mode squeezing) 

X (x) = exp {-2ixp) (3) 
in 



F^exp(^^(x2+p2)j (4) 

5(C) ^ exp (^^a^-^at^) (5) 

(C) = exp (C*aiaj - Cafa]) , (6) 

where a = (x + ip)/ \f2 and = {x — ip) /v^ are the ladder operators. Here the 
first three are single-mode operations and the last is our interaction. We describe 
these operations in more detail in section 4. For one-mode squeezing, Equation 
(5), the variable ( G C describes the orientation and amount of squeezing (de- 
scribed later) and on the annihilation operator a, S {() acts as 

5(C)^a5(C) =e+''fW+ie-''pW, (7) 

where x^^^ = x cos ^ -hp sin ^ and p^^^ = —x sin ^ cos 9 are the rotated quadra- 
tures. The last operation is the two- mode squeezing operator Si_j (C), which 
reduces the variance in the relative position and total momentum of the two 
modes being coupled. 



One way to easily visualise the effect of these operations on the phase space of 
the modes is using Husimi Q-function plots. For example, in Figure 1 we present 
the plots for the X (x) and Z (p) = exp {—2ixp) displacement operations. The 
Fourier transform, F in Figure 1, gives a rotation of 7r/2 in phase space. We see 
that we can easily generalise the Fourier transform to F (t) = exp (^it (a;^ + p^))- 
Figure 2 shows the effect of one-mode squeezing applied to the ground state. 
The variance in one axis is increased while the variance in the conjugate axis 
is decreased to compensate. Two-mode squeezing generates a similar effect, but 
between the pair-wise positions and momenta of the two modes. However, if 
we look at the relative positions and momenta instead, we see the correlations 
between the two modes. 
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Fig. 1. A Q-function plot of positive displacements of position and momentum from 
the ground state to other coherent states. The third Figure demonstrates the Fourier 
transform for a rotation of 7r/2. 




Fig. 2. One-mode squeezing apphed to the ground state. Shown is how r and 6 param- 
eterise the squeezing in phase space. 



To achieve general universal CVQC, we need to go beyond the realms of 
Gaussian states (which are the eigenstates of the class of Hamiltonians only 



quadratic in x and p). To do this requires a non-Gaussian operation, such as a 
non-linearity (see the above mentioned Kerr Hamiltonian) , a cubic-phase gate 
[16], or simply measurement - although without a clever scheme this may not be 
very useful for computation due to the probabilistic nature of the process. 

Having established a set of universal gates for computation, we briefly de- 
scribe how to carry out a computation. Wc can create any eigenstate of an 
operator which can be generated by our universal set (equations 3-6 plus a non- 
linear operation). To do this wc first initialise to an appropriate coherent state 
(say the ground state |0)) and perform the appropriate sequence of operations 
which generates the desired unitary operation. We have shown (equations 1 & 
2) that we can generate any desired unitary operation from the space of those 
available to us starting with our basis set but in practice there could be a long 
sequence of elementary operations to achieve the desired operation. After apply- 
ing the desired operation we can use homodync-likc measurements to find the 
output state and hence determine the result of the computation. The feasibility 
of an individual computation is at the mercy of the comparison of the resources 
available and the resources required, such as coherence times, number of modes 
etc. 



3 What is a micromaser? 

To perform CV quantum computations we propose the use of a micromaser-like 
system. The micromaser [17], or microscopically pumped maser, is essentially 
a very high quality microwave cavity with a rarified beam of atoms passing 
through it. The transition of single atoms moving through the cavity one at a 
time from one highly excited Rydberg state to another will keep the microwave 
field in the cavity primped. This is an extreme case of a beam maser, which would 
normally use a dense beam of atoms to pump a lossier cavity. The micromaser 
has been historically important as a test of cavity QED as a physical realisation 
of the Jaynes-Cummings model [18] which predicts the behaviour of a two level 
atom in a single mode field without semi-classical approximations. This model 
is important as the foundation of our understanding of masers and lasers. A 
simplified schematic can be found in Figure 3. 

By virtue of the very high quality of a micromaser cavity, for which the best 
examples may retain a photon for 0.3 s [19], the linewidth of the cavity mode is 
very small. For cavity of lifetime 0.3 s the linewidth is less than a hertz for a mode 
frequency of 21.456 GHz. The field in a micromaser cavity is well suited to CVQC 
due to this long lifetime and well defined mode frequency. Quantum coherence of 
the cavity field is also maintained by the high quality cavity. We propose the use 
of the quadratures of the cavity field as a continuous variable for CVQC. Multi- 
variable computations require extensions of the micromaser to multi-mode fields 
which requires modified cavity designs and carefully chosen energy levels of the 
atom. The Rydberg atoms that are commonly used in micromaser experiments 
provide a multitude of possible transitions to couple to many modes. 




Fig. 3. The 'phase sensitive' micromaser. Atoms are emitted by the source in the 
excited state and may be rotated coherently between the ground and excited states 
(both of which are actually highly excited Rydberg states of the atom) with a lossy 
microwave field, labelled as rotation here. The cavity in the centre may then interact 
with the atom, and a final rotation on the atom allows a measurement basis to be 
selected for measurement of the atomic state using state selective field ionisation. 



3.1 The Jaynes-Cummings Model 

The Jaynes-Cummings model describes the simplest non-trivial interaction be- 
tween atoms and light, consisting of a single atom with two states, labelled |e) 
for the upper and \g) for the lower^, interacting with a single mode field. The 
Hamiltonian of this interaction is [20] 

H = —^(73 + hioa^a — ihg (^a~^d — &~a') , (8) 

atom field interaction 

where fkUa is the transition energy of the atom from \g) to |e), w is the frequency 
of the mode, g is the atom- field coupling constant, a^{a) is the creation (annihi- 
lation) operator of the field mode and (T+((T~) is the atomic raising (lowering) 
operator. This model can be generalised easily to include more field modes or 
atomic levels. 



3.2 Some experimental considerations 

The source of atoms used in a micromaser is typically an effusive oven. The 
atoms emitted will thus have Poisson distributed arrival statistics and a velocity 
distribution determined by the oven temperature. The velocity distribution can 
be refined by the use of Fizeau wheel or detuning the first step excitation laser 
and placing it at an angle to the atomic beam [21] to select the velocity which 
Doppler shifts the laser into resonance. Poisson statistics of the atomic beam 
mean that there is some chance of multiple atom events in the cavity which are 
undesirable. The strategy used is to rarify the beam such that the chance of 
two or more atom events is lowered. The probability of a one-atom (desirable) 

These stand for excited and ground respectively, although in a laboratory two highly 
excited Rydberg states of rubidium are chosen such that the difference in energy 
between the states is close to resonance with the microwave cavity field. 



event is given by Pi = q-'^'^^I^ [21], where r is the rate of atomic arrivals at the 
cavity, L is the length of the cavity and v is the velocity of the atoms. Given 
that a typical micromaser cavity is approximately 3 cm long and that the selected 
velocity of the beam will be approximately 300 m s~^ , an average rate of 10 atoms 
per second will result in the probability of an event being a one atom event of 
99.8%. Operations that require many atoms may have to sacrifice some precision 
due to the the greater atomic rates necessary to complete a computation within 
the decoherence time of the cavity field. For experiments run by the micromaser 
group in Garching, cavities as good as up to a quality factor of Q = 4 x 10^° 
were used, corresponding to a life time of the field of 0.3 s. 

Some of the most recent advances in micromaser physics have involvcid the 
generation of Fock (number) states [22]. Using a micromaser field in a Fock state 
it is possible to produce a single atom in a particular state on demand [19]. If 
atoms on demand can be supplied then the rate can be greatly increased as 
multiple atom events are effectively eliminated, and many more computational 
operations performed. Another intriguing possibility for a single atom source is 
with the use of a standing-wave dipole trap, which can be used to accelerate 
single atoms deterministically [23]. 

4 Operations in the Micromaser 

Having described the two main schemes for encoding CVs in quantum systems for 
computation in Section 2, we need to choose a scheme for use in the micromaser. 
To decide this, it is necessary to look first at which states we can initialise 
and what operations we can perform on them. Goherent states^ [20], number 
(Fock) states (including Trapping States) [22] , steady states [24] , and tangent & 
cotangent states [25] may all be produced in a micromaser cavity field. Most of 
these are unsuitable for UCVQG but we see that we can create (;oherent states. 
We now show how we can achieve UCVQG on the Gaussian states, going through 
the list of required operations. 

4.1 Displacement Operations 

The single-quadrature displacement functions 

X {x) = (T^'^'P , Z{p) = e^^P* (9) 

are of the form of the generalised Displacement Operator, 

D{a = x + ip) = exp {2ipx — 2ixp) . (10) 

Explicitly, X {x) = D [x) and Z {p) = D (ip). These operations are depicted in 
Figure 1. 

® We use the standard definition of a coherent state, \a) = e~'"' ''^ JDn I"')- 



Displacement is remarkably simple to achieve in the micromaser cavity field 
[26, 27]. By applying the appropriate external coherent field we can displace the 
state in the cavity by D (a), thus by choosing an appropriate a, we can achieve 
X {x) and Z (p) in the micromaser. The experiments of Lange and Walther 
[27] demonstrated control over the average photon number over a large domain, 
from the subphoton level for states less than a = 100. Via feedback we can make 
these as stable as we require [28]. With modern microwave synthesisers much 
better performance is expected, however without knowledge of the device used 
it is difficult to estimate this improvement over the performance in the original 
experiment. In an actual experiment the power required for specific shifts should 
be determined experimentally as each cavity will couple at a diff'erent strength 
to an external field due to differences in machining. 

4.2 The Fourier Transform 

The Fourier transform F = e''^(*'+P^)/2 = ^^^(^+1)/^ jg simply a 7r/2 rotation 
in phase space about the origin. For example, on a ground state displaced in 
the positive-i; direction, we obtain a ground state displaced in the positive-p 
direction, as shown in Figure 1 in Section 2. It has a trivial class of eigenstates, 
the number states: 

F \n) = e''^("+i)/2 |n) = jn) . (11) 

This can easily be achieved using a micromaser with the atom detuned from 
the mode which we wish to apply the Fourier transform to, and far detuned 
from the other modes so as not to act on those. Given that the linewidth of each 
mode is so relatively small, this is a trivial requirement. For a mode detuned by 
A from an atomic transition the atom does not make a transition from |e) to 
1^), but the combined system evolves for a time t as 

|e,n) ^e-'»'("+i)*/^|e,n) , (12) 

and the atom may be neglected after interaction. The detuning can clearly be 
chosen to satisfy the Fourier transform in Equation (11). 

4.3 One-Mode Squeezing 

As there is no direct squeezing operation in the micromaser, we will first describe 
squeezing for linear optics since the formalism is the same. In linear optics, one- 
mode squeezing 

5(C)=exp(^a2-^at2) (13) 

is generated via a nonlinear-optical x'^^ interaction and yields the attenuation 
of a quadrature and the amplification of its conjugate. This is called nonlinear 
in optics, but to us this is a linear operation since it's only quadratic in the field 
operators. We'll still need a higher-order nonlinearity for our universal gate set. 



Since we are not considering linear optics, an analogue must be found for a 
micromaser system. The simplest system to consider is the two mode micromaser. 
This is best described as a three level ladder of atomic states for which the 
difference in energy between the uppermost and the lowermost states is twice 
the frequency of the field. The middle state is detimed from the one photon 
transition. As the central state is detuned from resonance with the field, the 
effective process is two photon transitions between the upper and lower state of 
the atom. The atom must be prepared in a particular superposition and after 
interaction the atom must be pulsed with a classical field so that it does not 
'give away' information about the cavity field and lead to decoherence. 

The particulars of this system were discussed in a paper by Orszag et. al. [29], 
but it is intuitive that Hi = k (a^ — a'l'^) will be the form of the two photon 
process. By analysing the cavity's steady state, if it begins in a superposition 
^„A„|n), then it tends towards a squeezed state if \a/j\ — 1 — e for a very 
small e. This state is highly squeezed in momentum but is not quite of minimum 
uncertainty. We cannot take an infinite limit as the amount of squeezing is 
dependent on e and all the squeezing disappears for e = but we can got an 
arbitrary amount of squeezing. After many atoms have passed through the cavity, 
we approach a pure squeezed state with ~ 0. 

Since it can take a long time to achieve this amount of squeezing, we may wish 
to prepare a sufficient quantity of squeezed states in advance in order that we 
may transfer the squeezing onto a state during the computation. Lam et. al. [30] 
have demonstrated this using Optical Parametric Oscillators in linear optics 
with the crucial component of a beam splitter to couple the modes. We show in 
Section 4.4 that a beam splitter may be replaced by any interaction and then 
explain an achievable interaction for the micromaser. 



4.4 Two-Mode Squeezing 

In linear optics, two-mode Squeezing 

Sij (C) = exp (C*aiaj - Caja]) (14) 

is generated via a non-degenerate optical parametric amplifier (NOPA). It gen- 
erates correlations between two modes, so it is an entangling operation. 

While it is a simple operation in mathematical terms, it is not as easy to 
achieve in the micromaser. However, the result for universality given in Section 
2 tells us that we need any interaction between two modes. By passing a 3-level 
atom with the energy level diagram in Figure 4 through a cavity with two modes 
we can get an interaction between the modes which effectively excludes the atom. 
The interaction also involves a coherently pumped lossy field and allows photon 
transfer between the two relevant modes through virtual processes. It is governed 
by the effective Hamiltonian with the assumption that 61,62,63 ^ gi, g2, T 



Hes = 0+ ^^p^ (a\al - 0162) , 



(15) 



i jf_ 

\a) 

Fig. 4. The energy level diagram for the 3-level atom. (5„ is the detuning of the relevant 
mode from each transition, and cOn is the frequency of each mode. The thick arrow 
between \a) and |c) denotes a classical pump field. If the atom moves up the ladder 
from |o) \b) |c) it will then be pumped back down to \a) by the classical field, 
acting on the field as the net operator 01O2. The converse is also true so that this 
system will behave like two-mode squeezing. 



where acts as a Fourier transform component on each mode. The rest is clearly 
related to two mode squeezing as the Schrodinger equation is solved to give the 
evolution operator U = c"'^"^'*. The coupling strength gi (92) is between the 
first (second) cavity mode and the transition between states |a) and \b) (|6) and 
|c)) of the atom, F is the coupling strength between the coherently pumped lossy 
field and the transition between levels \a) and |c) and 

p=J2alai(^^^d?j, Q = P + 62. (16) 

In an actual experiment this final transition would likely be replaced with a two 
photon transition to follow selection rules. We need to modify the state vec- 
tor simultaneously by j-!/)') = exp (^i'^^^i 2 (wj + 6i — ^) Now, two 
modes inside a single cavity can be coupled and we have an interaction polyno- 
mial which is good enough.. The effective coupling constant of the interaction is 
9ig2r/PQ and the interaction is very close to two mode squeezing. 

4.5 Non-linearity 

As well as the simple linear and interaction terms listed above, we need some 
non-linearity to achieve universal (not just Gaussian state) computation. We can 
simply measure the state of a 2-level atom (see Section 3) after transit through 
the cavity. This does not give a clean Hamiltonian, but any non-linearity is 
needed to be able to generate any Hermitian polynomial of continuous variables. 

5 Summary and Future Work 

In this paper we have given an overview of the history of continuous variables 
in computation, both classical and quantum. We stated the standard result for 



universal continuous variable quantum computation (CVQC) for two different 
encodings. We then gave an account of the micromaser, both what we can do 
in the experiment and what states can be created in the cavity. Given this, we 
showed how we can in principle achieve universal CVQC using Gaussian states 
in the micromaser using simple interactions. 

Having a system which can, in principle, perform universal CVQC is very 
useful but we must consider how feasible it is to use the micromaser in such 
a way. All the necessary operations can be produced quickly and accurately in 
the system and any Hamiltonian can be generated by a polynomial number of 
the base operations. We may have a problem of scalability since our current 
interaction is between two modes in one cavity. However, for low numbers of 
modes the micromaser is a perfect candidate for efficient universal CVQC. 
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